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We investigate the chemical composition and adhesion of chemical vapour deposited thin-film 
alumina on TiC using and extending a recently proposed nonequilibrium method of ab initio ther- 
modynamics of deposition growth (AIT-DG) [Rohrer J and Hyldgaard P 2010 Phys. Rev. B 82 
045415]. A previous study of this system [Rohrer J, Ruberto C and Hyldgaard P 2010 J. Phys.: 
Condens. Matter 22 015004] found that use of equilibrium thermodynamics leads to predictions of 
a non-binding TiC/alumina interface, despite the industrial use as a wear- resist ant coating. This 
discrepancy between equilibrium theory and experiment is resolved by the AIT-DG method which 
predicts interfaces with strong adhesion. The AIT-DG method combines density functional theory 
calculations, rate-equation modelling of the pressure evolution of the deposition environment and 
thermochemical data. The AIT-DG method was previously used to predict prevalent terminations 
of growing or as-deposited surfaces of binary materials. Here we extent the method to predict sur- 
face and interface compositions of growing or as-deposited thin films on a substrate and find that 
inclusion of the nonequilibrium deposition environment has important implications for the nature 
of buried interfaces. 



I. INTRODUCTION 

Interfaces and surfaces are present in practically all de- 
vices and their detailed structure is typically crucial for 
the overall device functionality p~H6] . Understanding and 
ultimately controlling the thin-film deposition, the chem- 
ical composition and the adhesion at interfaces [7HT8] is 
a very desirable goal of industrial and scientific research. 
Characterisation of atomic structure and binding at in- 
terfaces is a fundamental step towards this goal. 

Atomistic modelling of materials [19-27 using meth- 
ods based on ab-initio density functional theory (DFT) 
allows for a detailed understanding of structure at the 
atomic and electronic level. Developing reliable mod- 
elling methods is of particular value for characterisation 
of interfaces which are buried inside materials and there- 
fore difficult to characterise experimentally with atomic 
resolution [20[ I28H5T] . Atomistic modelling has a poten- 
tial to accelerate innovation, for example, in the develop- 
ment of coatings and in the design of functional surface 
and interface materials [32]. A key element of a reliable 
method is a proper treatment of thermodynamic effects of 
a surrounding environment during creation of thin films 
and interfaces. 

Until recently, ab-initio thermodynamics (AIT) meth- 
ods were essentially methods of surface equilibrium (and 
here denoted as AIT-SE). These methods have focused on 
oxide surfaces [33j [34] or metal/oxide interfaces [35] as- 
suming equilibrium between the oxide surface (interface) 
and O2 in an 02-dominated, e.g., ambient environment. 
Oxide surfaces are in direct contact with the environment 
and the oxygen content of this environment will therefore 
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have a strong influence on the termination of the oxide. 
However, it is by no means clear how a surrounding could 
easily influence the composition at interfaces (which are 
by definition buried and insulated from the gaseous en- 
vironment). In fact, we have shown [36] that adapting 
and applying this AIT-SE method to the TiC/alumina 
interface predicts a structure and composition that pos- 
sesses essentially no binding across the interface. This 
result is evidently in conflict with the actual use of chem- 
ical vapour deposited (CVD) TiC/alumina multilayers as 
wear-resistant coating on cemented-carbide cutting-tools 

[to Eg. 

The present paper demonstrates that realistic descrip- 
tions of deposition environments are crucial for charac- 
terising thin-film and interface compositions and, as a 
consequence, adhesion to the underlying substrate. We 
show how effects of a steady-state deposition environ- 
ment (for example CVD) can be embedded into atomistic 
modelling. We employ a method of ab-initio thermody- 
namics of deposition growth (AIT-DG) [39], but extend 
it here to suit the more complex problem of understand- 
ing thin-film formation and corresponding as-deposited 
interfaces. We focus our discussion on TiC/alumina. 

The fundamental strategy is to compare free energies 
of reaction associated with thin-film configurations that 
differ in their detailed chemical composition. The key el- 
ements of the method are: (i) analysis of free energies of 
reaction G r and (ii) use of rate-equation modelling for the 
pressure evolution of the deposition environment. The 
key variables that determine G r are the partial pressures 
of the various constituents of the environment. Assuming 
a steady state, all partial pressures can be expressed in 
terms of a few rates. No assumptions about equilibrium 
between some of the species enter into this analysis. We 
point out that a steady state does not necessarily corre- 
spond to a state, sometimes described as dynamic equi- 
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FIG. 1: (Colour online) Structure of interfaces between TiC 
and various thin-film alumina overlayers. Colour coding: Ti 
= black, C = gray, Al = blue (black, small balls) and = 
yellow (light gray). The overlayers can be sorted into classes 
of thin-films with different nature of interfaces (A, B, C) and 
possessing different deviations from the full AI2O3 stoichiom- 
etry (here, Al deficiency only). The interface class predicted 
by the equilibrium AIT-SE method is shown in the top pan- 
els. Its non-binding character is obvious from the geometric 
structure (left) and electronic density (right), and in conflict 
with the wear-resistance of TiC/alumina coatings. The set 
of bottom panels shows alternative classes of thin-film alu- 
mina on TiC. These show stronger adhesion at the interface 
and have correct nature according to industrial use. They are 
found unstable according AIT-SE but stable according to the 
present nonequilibrium account, AIT-DG. 



librium [40] , where the system is assumed to gain no free 
energy by the deposition of stoichiometric alumina. The 
method also describes the evolving system in a certain 
range where dynamic equilibrium is not maintained. 

The paper is organised as follows. In Section [IT] we 
summarise the results of a previous structure search and 
equilibrium-thermodynamic analysis for TiC / alumina in- 
terfaces [36]. This section highlights the inconsistency 
in adhesion properties of the predictions made with the 
AIT-SE method. In Section III we give a general mo- 



tivation for the use of Gibbs free energies of reaction 
to predict the prevalence of thin-film classes and chem- 
ical compositions at as-deposited interfaces. Section [IV 
presents a simple model that describes the environment 
used for chemical vapour deposition of alumina on TiC. 
Section [V] contains the details for the evaluation of Gibbs 
free energies of reaction for TiC/alumina interfaces. In 
Section |VT] we present and discuss our results. We sum- 
marise our work and conclude in Section IYIII 



II. BACKGROUND 

In Ref. [36j we have presented an ab initio struc- 
ture search for thin-film alumina on TiC. We have con- 
sidered thin-film alumina configurations with different 



thicknesses and various chemical compositions. The con- 
figurations can be characterised by two numbers t and 
Aai- Here, t is the film thickness in terms of the number 
of O layers and Aai = -ZVai — 2/3iVo is the number of 
excess (Aai > 0) or deficiency (Aai < 0) in the number 
of Al atoms. 

All configurations approach the full stoichiometry in 
the thick- film limit. They essentially differ in chemical 
composition only at the interface, surface, or both. We 
therefore, in general, group the configurations into in- 
terface classes according to their corresponding value of 
Aai- Interface class A corresponds to Aai = — 4, inter- 
face class B corresponds to Aai = — 2 and interface class 
C corresponds to Aai = 0. Other values of Aai were not 
considered. 

Figure [T] presents energetically optimised TiC/alumina 
configurations for all three interface classes. The set 
of top panels details of the atomic (left) and electronic 
structure (right) at the interface class A. Straightforward 
adaption and use of the equilibrium AIT-SE method [33]- 
135] identifies interface class A as thermodynamically sta- 
ble over a wide range of temperatures and O2 pressures, 
see Ref. [36] However, no appreciable adhesion is found 
at interfaces of this type. These interfaces separate into a 
TiC substrate covered with a full layer of O (TiC/O) and 
a fully stoichiometric alumina overlayer. The electronic 
structure shows that the electron density essentially van- 
ishes between the TiC/O and the alumina, ruling out a 
significant covalent binding. Furthermore, the Al ions at 
the interface relax into the first O layer above the TiC/O, 
ruling out significant ionic binding. 

In the present work we quantify the (lack of) adhe- 
sion by calculating the ideal work of adhesion as W a dh = 

(-^substrate H - -^alumina -^TiC/alumina) / Here -^substrate 

and ^alumina are the energies of the isolated relaxed sub- 
strate (TiC for interface class B and C or TiC/O for 
interface class A) and the isolated relaxed alumina film 
(with in-plane lattice parameters constraint to the sur- 
face lattice of TiC). ^TiC/alumina is the energy of the re- 
laxed interface and S is the area of the contact surface. 
For interface class A we find a nonquantifiable (vanish- 
ing) value of I^adh, comparable to the uncertainty in the 
force relaxation in the underlying GGA calculations (de- 
scribed in Ref. [36]) . 

The non-binding character at interface class A 
is clearly in conflict with the wear-resistance of 
TiC/alumina multilayers [37] . We attribute this conflict 
to the fact that alumina does neither maintain equilib- 
rium with O2 during deposition 1 nor can one expect that 
it reaches such equilibrium after being removed from the 
deposition chamber. This true for the present focus on 
thin films and even more so for thicker overlayers. 



1 At best the alumina may maintain dynamic equilibrium with 
a number of gases in the surrounding. However, as we will also 
discuss below, even this assumption is too optimistic for the case 
of CVD of alumina on TiC. 
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FIG. 2: Schematics of elementary growth processes in an 
ensemble of systems that can be organised into three differ- 
ent classes A, B and C. The members of these classes could 
be alumina films (of various thicknesses but) sorted by their 
balance of Al and atoms and therefore having different sur- 
faces and/or interfaces to the substrate. The resulting inter- 
face classes may exhibit a very different nature of adhesion. 
Combining chemical reaction theory with a rate-equation de- 
scription of the system allows us to predict the probabilities 
for realisation of class A, B and C during steady-state growth. 



According to the AIT-SE method, interface class B 
(see bottom left of Figure [T]) is predicted to be stable 
only under extreme conditions, whereas interface class C 
(see bottom right panel) is predicted unstable over the 
whole range of allowed values of chemical potentials 
[36] . However, the works of adhesion associated with 
these interface classes, W a ah = 7.4 J/m 2 for interface 
class B and W a dh = 7.3 J/m 2 for interface class C, are 
in much better agreement with the wear-resistance of the 
material. We show below that the conditions prevailing 
during CVD of alumina allow for nucleation of either one 
of these binding interface classes B and C. 



III. AB INITIO THERMODYNAMICS OF 
DEPOSITION GROWTH 

In Ref. 39 we have introduced the AIT-DG method and 
demonstrated that the Gibbs free energy of reaction can 
be used as a predictor for the prevalence of the chemical 
composition at a growing surface. This was done for a 
binary material with two possible surface terminations. 
Here we extend this nonequilibrium description to the 
problem of identifying the composition of growing over- 
layers which exhibit both a surface and an interface to 
the substrate on which they are deposited. We follow 
a similar line of argumentation as in Ref. 39] and com- 
bine chemical reaction theory [40 , with a rate-equation 
description of the probabilities for finding any of the pos- 
sible film compositions. 

Figure [5] illustrates a collection of systems (e.g. an 
ensemble of growing thin films) that are grouped into 
three thin-film and interface classes A, B and C . The 



members of the classes can have different film thicknesses 
but are sorted according to their chemical composition in 
terms of excess or deficiency atoms of a specific species. 
Members of different classes are allowed to transform into 
members of other classes via the three reactions labeled 
as /, //, and all being characterised by forward and 
backward rates T^ and T l h (i = I, II, III). Growth of, 
for example, a film of class A results by a net flow along 
the reaction chain, A— ?>B— ?>C— ?>A. 

Chemical reaction theory relates the forward and back- 
ward rates to the inverse temperature j3 (in units of en- 
ergy) and the Gibbs free energy of reaction by (3AG l r = 
— lnTj/r^. At the same time, the probability for a ran- 
dom member of the ensemble to belong to one of the 
three classes is described in terms of rate equations, 
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This follows from a straightforward generalisation of the 
analysis presented in Ref. [39] The steady-state solutions 
for the probabilities can be expressed as ratios of sums 
of products of reaction rates, for example, 
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We use the differences in Gibbs free energies of reaction 
as an approximate predictor for the prevalence of the 
different classes of films and interfaces, 
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The evaluations of relative probabilities Q are exact in 
the limit where AG T r + AG 1 / + AG 1 / 1 = or, equiva- 
lently, T^/V 1 / 1 = r£r£ 7 r£ /7 . This limit corresponds 
to dynamic equilibrium in stoichiometric growth of thin 
films (irrespective to which interface class they belong). 
We note that these predictors [and the full evaluation 
([2|] are beyond a simple assumption of detailed balance 
because the probabilities are steady-state but not equi- 
librium probabilities. 

With any use of the presented AIT-DG modelling, one 
should always check the quality of the approximative pre- 
dictors Q . The appendix provides the theoretical frame- 
work for such a test in the case of alumina thin-film de- 
position and interface formation on TiC. 
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FIG. 3: Schematic model of chemical vapour deposition of 
alumina. A H2-AICI3-CO2 gas mixture with relative concen- 
trations ch 2 , CAICI3 and cco 2 is supplied at rate Rs to a hot 
reaction chamber which is kept at constant temperature. At 
the same time, there is a continuous exhaust of gases at rate 
Re, also keeping the pressure constant. Including the reac- 
tions (4a) and (4b ), taking place inside the chamber, we model 
this system by the set of rate equations The resulting 
steady-state partial pressures can be used to determine the 
free energies of reactions ^ for different TiC/alumina sys- 
tems. 



IV. MODELLING OF CVD OF ALUMINA 

A. Rate-equation modelling 

Figure [3] illustrates the CVD process utilised for alu- 
mina growth on TiC. A H2-AICI3-CO2 supply gas mix- 
ture with relative concentrations q is supplied to a hot 
chamber at rate Rs. Inside the chamber water and alu- 
mina form (in parallel) at rates Ru 2 o an d ^ai 2 o 3 accord- 
ing to [4T] 



H 2 + C0 2 ^> H 2 + CO, 



(4a) 

2AICI3 + 3H 2 ^ Al2 ° 3 ) A1 2 3 + 6HC1. (4b) 

The total pressure inside the chamber is kept at a con- 
stant value by a continuous exhaust at rate Re of both 
reaction products and unused supply gases. 

We describe the evolution of the environment by indi- 
vidual partial (ideal gas) pressures with a coupled set of 
rate equations, 



Pt 



d m a aR s - jRe + K 2U Ru 2 o + ^ a12U3 ^ai 2 o 3 . (5) 

Here, pi = Pi{t) is the momentary pressure of chemical 
species i inside the reaction chamber, p = p(t) = J2iPi(t) 
is the momentary total pressure, q is the concentration 
of the chemical species i in the supply gas, and i^ 2 ° and 
V f 2 ° 3 are the stoichiometric coefficients 2 of the chemical 
species i in reaction Hal) and (4b), respectively. 



The resulting steady-state partial pressures (dtPi(t) = 
0? p(t) = p = const.) can be expressed in terms of 



scaled reaction rates rn 2 o = Rn 2 o/Rs and rAi 2 o 3 = 
Rai 2 o 3 /Rs, 

Pi = P f- — • (6) 

1 + ^Al 2 O s 

This result allows us to determine free energies of reac- 
tion (see below) from a few experimentally controllable 
(p, T, Rs) or at least measurable (Rn 2 o, Ra\ 2 o 3 : Re) 
quantities. 



B. Limits on reaction rates 

The reaction rates possess natural bounds that cannot 
be exceeded in steady state. First, for each reaction an 
upper bound for the reaction rate is given by the con- 
dition of dynamic equilibrium in this reaction. Reaction 
(4a) can proceed from left to right only if the chemical 
potentials fulfil /iH 2 + Mco 2 > Mh 2 o + Mco- Similarly, re- 
action (4b) requires that 2/zaici 3 + 3/iH 2 o > AGai 2 o 3 + 



6/iHCb where AGai 2 o 3 is the gain in free energy per sto- 
ichiometric formula of Al 2 03. Dynamic equilibrium in 
these reactions is reached if equality holds which corre- 
sponds to a specific values of Ru 2 o and Ra\ 2 o 3 - 

In the present case another bound is found by compar- 
ing reactions (4a) and (4b). The latter reaction requires 



three units of H 2 0, while the former produces only one. 
Thus, in steady state, Ra\ 2 o 3 < 3i?H 2 o = ^ai 2 o 3 - How- 
ever, we do not expect our model to be applicable in 
the the limit where Rai 2 o 3 ~^ ^ai 2 o 3 - ^ n this am ~ 
mina deposition becomes instantaneous and on average, 
the water pressure will vanish. The ideal-gas descrip- 
tion of the environment will therefore be inappropriate 
and kinetic aspects become dominating. We choose to 
consider situations where Ra\ 2 o 3 is sufficiently separated 
from it— 3 - 



V. AB INITIO EVALUATION OF GIBBS FREE 
ENERGIES 

The Gibbs free energies of reaction are calculated as 
follows. We consider a general AImOat film as the prod- 
uct of a (hypothetical) chemical reaction starting from 
the substrate (where M = N = 0). Reaction ( [4b] ) de- 
scribes stoichiometric solidification of alumina. In addi- 
tion, the CVD environment enables the following reaction 
pathways for non-stoichiometric deposition, 3 



A1C1, 



3 
2 

H 2 



H 2 -> AW + 3HC1 (7a) 

Ofilm+H 2 . (7b) 



2 Stoichiometric coefficients are counted negative if a species is 
consumed and positive if a species is produced in a reaction. 



We do not include the reaction CO2 — > O eX c + CO since the 
associated free energy of reaction is higher than that associated 
with ffb). 
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The label 'film' indicates the incorporation into the alu- 
mina film on the substrate. 

We define the free energy of reaction, G^ ,N , associ- 
ated with the (hypothetical direct) deposition of a gen- 
eral AImOat film, 



G 



M.N 



Gt\c/a\ m o n - Gxic 



iV(/iH 2 -/iH 2 o)- 



MAlCl f 



(8) 



Here, G T ic/Ai M o N 1S the free energy of the film adsorbed 
on the TiC substrate, GxiC is the free energy of the clean 
(relaxed) TiC substrate. The set of chemical potentials 
[ii describe the Gibbs free energy of the various gases 
that contribute to the reaction. 

The Gibbs free energy definition ([8| is unambiguous 
and relevant for our nonequilibrium thermodynamic ac- 
count of thin-film growth and interface formation. Of 
course, no direct reaction for the deposition described by 
J8|) exists. However, stoichiometric combination of (7a) 



and (7b) effectively reduces to (4b). We treat G^'^as 
a conservative function of the state variables M and N 
(that is independent of the details of the order of the 
deposition steps that result into an AImOat film). As a 
consequence, we can extract the values of AG r associated 
with the reaction that transforms a film characterised by 
M r and N r into a film that is characterised by M and 
N as difference between G?> N and Gf > N ' . These dif- 
ferences allow the direct calculation of the relative prob- 
abilities in the limit of dynamic equilibrium from our set 
of predictors (J3|. 

We consider alumina films of the type Al4 t _406t (inter- 
face class A), Ai4t-206£ (interface class B) and Al^Oet 
(interface class C), where t is the thickness in terms of 
the number of O layers. Two times the reaction (7a) thus 
corresponds to reaction / and in Figure [2] and six times 
the reaction (7b) corresponds to We note that, al- 
though / and II are both formally described by the gas 



reaction (7a), they are different and possess different free 



energies of reaction since the substrates on which they 
occur and the final products are different. 

We evaluate Gibbs free energies of reaction G^' N as 
described in Refs. |33] [34j |39j The free energies of solid 
phases (substrate and potential thin films) are replaced 
by their DFT total energies, G S oiid ~ ^soiid- 4 For gases, 
we employ the ideal-gas approximation, 

M(T, Pi ) = a + A°(T) + k B T\n( Pi / P °). (9) 



4 In fact, we correct the total energies of the films by subtract- 
ing the strain energy of the stoichiometric part of the film, 
Em m -> Emm ~ n A i 2 o 3 A strain , where A strain is the is differ- 
ence between the strained (due to the expansion to the TiC lat- 
tice in the interface plane) and the unstrained bulk alumina per 
stoichiometric unit. 



0> 

o 

> 




-6 -5 -4 

10 SlO r H.O 



FIG. 4: Gibbs free energies of reaction G r for formation 
of CVD thin-film alumina with different thickness and dis- 
playing different interfaces classes as functions of the scaled 
reaction rate rn 2 o- The vertical line corresponds to dynamic 
equilibrium in the water formation reaction. Parameters of 
Ref. [41] have been used for the composition of the environ- 
ment. A deposition temperature of T = 1000 °C, a pressure 
of p = 500 mbar and rAi 2 o 3 = th 2 o/3.1 were assumed. No 
qualitative changes arise when these parameters are varied in 
range typical for CVD of alumina [JT]. The films are labeled 
according to their corresponding interface class (A, B or C) 
and their thickness in terms of O layers (here the cases for 3 
and 4 layers are displayed). 



Here e$ is the DFT total energy of the gas (molecule) and 
A^(T) is the temperature dependence of /i^ at a fixed 
pressure p°. We use the values tabulated in Ref. [43] for 
p° = 1 atm. For the individual partial pressures pi we 
use steady-state pressures (loT) specified by CVD process. 



VI. RESULTS AND DISCUSSIONS 
A. Gibbs free energies of reaction and growth 

In Figure [4] we plot the Gibbs free energies of reaction 
G r [see ^] for various thin- film alumina configurations 
on TiC. The films are labeled according to their corre- 
sponding interface class (A, B, C) and their thickness. 
We assume typical values for the the deposition temper- 
ature T, the total deposition pressure p and for the con- 
centrations Ci of the different gases in the supply gas [41] : 
T = 1000°C, p = 500 mbar, c A ici 3 = 0.04, c C o 2 = 0.04, 
and H 2 constitutes the balance. 

The vertical line corresponds to dynamic equilibrium 
in the water-producing step (4a). We emphasise that, in 
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general, dynamic equilibrium in the water-producing step 
does not imply dynamic equilibrium in the alumina de- 
position. In the figure we have chosen rAi 2 o 3 = r H 2 o/3.1. 
For this value of rAi 2 o 3 the dynamic equilibria roughly 
coincide within our our approximation for the Gibbs free 
energy variation [36]. This follows from the observation 
that the values of G r for the B3 and B4 films (which 
differ in thickness by one full layer or two stoichiometric 
units) are approximately equal on the vertical line. We 
note, however, that there is an uncertainty in this value 
°f r Ai 2 o 3 f° r coinciding dynamic equilibria. The reason is 
the uncertainty in the calculation of chemical potentials 
(in particular total energies of the molecules, see discus- 
sion in Ref. [39]) and the uncertainty in the total energies 
of the films themselves. For the latter to be accurate we 
would have to make sure that their atomic structures are 
fully optimised. The geometries identified in Ref. [36] and 
used here are candidates for the optimised structures but 
not guaranteed to the optimised structures. 

Figure 4] shows that, within the possible range of rH 2 o, 
all films have a negative value of G r . As a consequence 
all films are stable with respect to the substrate. This 
remains true also if we decrease the value of rAi 2 o 3 - 

The figure also demonstrates stoichiometric growth of 
alumina films. Stoichiometric growth corresponds to 
AG* + AG 1 / + AG 1 / 1 < 0. This condition is fulfilled 
as can be seen from the fact that the B4 film has lower 
free energy of reaction than the B3 film and that the for- 
mer consists of two stoichiometric units of alumina more 
than the latter. 

We emphasise that the Gibbs free energy variations 
shown in figure [4] directly reflect the closely related na- 
ture of process / and II. Process / and II differ only 
in the solid reactants and solid products, the gaseous 
reactants and products are the same. In chemical reac- 
tion theory [40] we can express the Gibbs free energy of 
reaction as (3AG r = lnT^/r^ = — lnif eq + InQ. Here 
K eq = fcf/fcb is the equilibrium constant of the reaction, 
fcf (fcb) is the forward (backward) reaction rate constant 
and Q is the reaction quotient (ratio of concentrations 
of gaseous products and of reactants). Since the gaseous 
reactants and products in process / and II are identical, 
the reaction quotients for reaction / and II are identi- 
cal. Therefore, /3(AG r J - AG 1 /) = ln(K^/K^) must be 
constant. In figure [4] we have 



B. Thermodynamic analysis 

Figure [5] reports the calculated predictor ^ for the 
prevalence of interface class A and C relative to inter- 
face class B as functions of the scaled reaction rate for 



water formation [see (4a)]. This predictor corresponds to 
the approximate relative steady-state probabilities (being 
exact in the limit of dynamic equilibrium in the alumina 
deposition). We have tested the quality of our predictor 
over a wide range of possible choices for unknown rate 
constants; we refer to the appendix for a more detailed 
presentation. 

The vertical line corresponds to dynamic equilibrium 
in the water-producing step (4a). We assume the same 
values as before for the temperature and the concen- 
trations of the different species in the supply gas [4T] . 
For the total deposition pressure and the reaction rate 
of alumina deposition we consider the following pairs 
of parameters, p = 500 mbar and rAi 2 o 3 = ^h 2 o/3.1, 
p = 10 mbar and rAi 2 o 3 = ^h 2 o/3.1 and p = 500 mbar 
and rAi 2 o 3 = nH 2 o/300. The location of the vertical line 
is independent of these parameters (but depends on tem- 
perature). 

We find that, in or close to dynamic equilibrium in 
the water producing step, the (approximate) probabil- 
ities for the prevalence of interface class A and C are 
much smaller than the probability for the prevalence of 
interface class B, Px << Pb (X = A or C). As rH 2 o de- 
creases Pq increases and becomes larger than Pb at some 
point. The value of Pa, on the other hand, decreases as 
rn 2 o decreases. Thus, our results show that the nonbind- 
ing interface class A is not realised in CVD of alumina 
as described in Section IV Instead, if reaction (4a) takes 



place sufficiently close to dynamic equilibrium, interface 
class B is predicted. 

We note that it is by no means obvious that dynamic 
equilibrium is always maintained. If that was the case, 
r Ai 2 o 3 would always assume its maximum value. Then, 
however, it would be possible to scale the absolute de- 
position rate Ra\ 2 o 3 = r Ai 2 o 3 * to infinity simply by 
increasing the supply rate. It is against common sense to 
expect the deposition to remain in dynamic equilibrium 
as we increase the supply flux. 



Characterisation of interface adhesion 



P(AG{ - AG 1 /) 

= [G r (P4) - G r (A4)] - [G r (G3) - G r (P3)] 
= [G r (P4) - G r (G3)] - [G r (A4) - G r (B3)]. 



(10) 



We notice that the curves corresponding to B4 and C3 
films are approximately parallel. The same applies for 
the curves corresponding to A4 and B3 films. Thus, both 
differences after the second equal sign in ( 10 ) are constant 



and the AIT-DG results can also be used to compute 



Figure [6] details the electronic density at the interface 
class B. We have chosen the same isosurface level as in the 
top right panel of figure [I] The Ti-0 bonds show both 
ionic and covalent character. The density also reveals the 
partial covalent character of the otherwise highly ionic 
alumina itself. This covalent character can be seen by a 
small amount of electrons spilling over from neighbouring 
O layers. 

A Bader analysis including core charges [44U46] shows 
that the ionicity of the Ti atoms is lower at the interface 
B (#Ti = 1-71 e) than at the interface A (gxi = 1-88 e, 
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top panels of figure [TJ. At the same time, the ionicity 
of the O atoms is higher (qo = —1.46 e) at interface 
B than at interface A (q — —1.13 e). The Ti-0 bond 
strength is lower at the interface B than at the interface 
A. This is reflected by an increased Ti-0 layer separation 
at interface B (/xi-o = 1.48) with respect to that at in- 
terface A (7xi-o = 0.88). However, the increased ionicity 
of the first O layer at interface B and the presence of Al 
ions between the first and the second O layers above the 
interface plane give rise to a strong overall interface ad- 
hesion. This adhesion is absent at interface A. Thus, the 
strong binding at the wear-resistant interface B that is 
most likely to be created in CVD of alumina arises from 
an interplay between a softening of the Ti-0 bonding, an 
increased ionicity in the first O layer and a strong ionic 
binding in the alumina. The latter is supplemented by 
weaker covalent bonds between O layers in the alumina 
coating. 



FIG. 5: Predictors for prevalence of interface class A (solid 
lines) and interface class C (dashed lines) as a functions of the 
scaled reaction rate ru 2 o for water formation. The predictors 
are the approximate logarithm of the probability relative to 
the probability of prevalence of interface class B. The verti- 
cal line limits th 2 o to the right and corr espo nds to dynamic 
equilibrium in the water formation, see (4a). Parameters of 



Ref. ^2 have been used for the composition of the environ- 
ment and a deposition temperature of T = 1000 °C was as- 
sumed. Black thick lines correspond to a deposition pressure 
of p — 500 mbar and rAi 2 o 3 = th 2 o/3.1, thin light lines to 
p = 10 mbar and rAi 2 o 3 = th 2 o/3.1, and thin dark lines to 
p = 500 mbar and rAi 2 o 3 = th 2 o/300. Close to dynamic equi- 



librium limit in reaction (4a), we predict the highest preva- 
lence for interface class B. As rH 2 o decreases, the likelihood 
for the occurrence of interface class C increases. This like- 
lihood strongly increases further with decreasing rAi 2 o 3 and 
also as the total deposition pressure increases. The nonbind- 
ing interface class A is highly unlikely to be created in the 
CVD process. 




FIG. 6: (Colour online) Electron density at the interface B 
which shows the highest probability to be created in CVD of 
alumina on TiC. Same colour coding as in Fig. [I] 



VII. SUMMARY AND CONCLUSION 

We have extended a recently proposed nonequilib- 
rium method for ab initio thermodynamics of deposition 
growth from application to surfaces terminations [39] to 
prediction of thin-film and interface formation. Our re- 
sults demonstrate that a careful treatment of deposition 
conditions in thin-film and interface modelling is crucial 
for understanding adhesion properties for CVD alumina 
on TiC. Assuming equilibrium between the oxide and 
O2 (as in [33-35 ) results in prediction of a configura- 
tion that shows no binding across the interface (interface 
class A) [36 , see top panel in Fig. [I] This is in con- 
flict with the wear-resistant nature and industrial use of 
the material [37] . In agreement with the wear-resistance 
of TiC/alumina coatings, we predict the deposition of 
strongly binding interface of type B (in or close dynamic 
equilibrium in the water forming step) or interface of type 
C (away from dynamic equilibrium). 

We expect that a similar analysis will be necessary also 
for other buried interfaces that form during a deposi- 
tion process in an environment that strongly differs from 
ambient conditions. The thermodynamic method used 
here only makes reference to the molecular species that 
are present (and directly relevant) during deposition. In 
contrast to equilibrium-thermodynamics approaches [33]- 
[35], this method is therefore not limited to oxides (al- 
though we here illustrate the method for a particular ox- 
ide) or materials that contain a constituent X for which 
a dimer X2 could serve as reference. Furthermore, the 
method allows for a search of conditions (supply gas com- 
positions, deposition temperatures, deposition tempera- 
tures) that favour deposition of a pre-specified interface 
composition with desirable properties. Our results sug- 
gest that the nonequilibrium ab initio thermodynamics 
method (Ref. [39 and present extension) can be useful 
in guiding experimental optimisation of present-day ma- 
terials and design of novel such. 
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Appendix A: Determination of film prevalence 

The presentation of our results for the prevalence of the 
different alumina films rests on the analysis of the predic- 
tors Q . These predictors are strictly valid only under the 
assumption of dynamic equilibrium (in the alumina de- 
position) . In this appendix we give a prescription of how 
the quality of these predictors can be tested with respect 
to the exact steady-state probability for the prevalence 
of different film classes. 

We consider the exact steady-state ratio between the 
probability for A and B, see equation We intro- 

duce Xi = Tl/T) - — ' prn " — 1 ~ - --i-- 1 -- 1 ™- 



x 3 \ The 



exp(/3G*) and 

functions a and X{ depend on rn 2 o but are otherwise 
completely determined by our ab initio calculations (and 
the thermochemical data to determine the chemical po- 
tentials). In particular, the parameter a measures the 
departure from dynamic equilibrium; a = 1 corresponds 
to dynamic equilibrium, a > 1 corresponds to growth 
beyond dynamic equilibrium and a < 1 corresponds to 
evaporation. We furthermore introduce F\ = ^f^ 1 / an d 
F 2 = r^r^. These functions also depend on rn 2 o but 
require some additional parametrisation, see below. 

Independent of the par ametrisat ions, we can express 
the ratio of probabilities Q as 



The quality of the appr oxim ate p redictors (|3| can be 
teste d by eval uating (Al) and (A2). The evaluation of 
( |A1[ ) and ( |A2[ ) requires the specification of the functions 
F\ and F 2 . Here, we assume the following parametrisa- 
tions, 



*1 = 



F 2 



1 / 
1 f 

1 / 

r 7 
1 / 



UII 



^\x{ u \ 



k f 



(A3) 



(A4) 



where the v\ can but do not need to be identical with the 
stoichiometric coefficients. We note that this parametri- 
sation in only a crude approximation and that growth 
can, in general, be more complicated. 

For alumina growth, the reactions / and II are for- 
mally identical. The only difference is that they take 
place on different substrates and generate different solid 
products. Therefore, F\ reduces to the ratio of the rate 
constants, which is constant. Reaction / and /// are 
different and F 2 therefore depends on rn 2 o- Thus we 
have 



Fx 



k 1 

UI 
K f 

UII 



[H 2 0r[H 2 ] 



k f 



[AlCl 2 ] r [H 2 ] s [HCl]- 



(A5) 
(A6) 



Pa 
Pb 



x x x 2 + x 1 F 1 F 2 + F 2 



-x 2 +F x F 2 



~ 1 x7 1 F 2 



(Al) where m, n, r, s, t are all positive and of the order of 1. 



For the ratio between the probability for C and B we 
find 



Pb 



ririn 
L b L b 



vyvi 11 + r I f r 1 / 



T n T in 



r / r // 

1 f L b 



FiF 2 + F 2 + axix 2 
x 2 F 2 + axix^ 1 + axix 2 FiF 2 



(A2) 



Since the rate constants kj, 



we have compared (3a) with (Al) and (3c 



k 1 / , and kj 11 are unknown, 



with (A2) 



for several choices of m, n, r, 8, t and a bro ad range 
of ratios of rate constants. We find that (A2) is essen- 
tially described by (3c). For (Al) the approximation 



(3a) can deviate as rH 2 o decreases but the qualitative 
result (Pa << Pb) is n °t affected. Moreover, the ratio 
Pa/(Pb + Pc) is always strongly suppressed. 
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